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Dimensional reduction is a key issue in finite temperature field theory. For example, when following the 
QCD Free Energy from low to high scales across the critical temperature, ultrasoft degrees of freedom can be 
captured by a 3d SU(3) pure gauge theory. For such a theory a complete perturbative matching requires four loop 
computations, which we undertook by means of Numerical Stochastic Perturbation Theory. We report on the 
computation of the pure gauge plaquette in 3d, and in particular on the extraction of the logarithmic divergence 
at order g 8 , which had already been computed in the continuum. 



INTRODUCTION 

Numerical Stochastic Perturbation Theory 
(NSPT) methods .1 arise within the Stochas- 
tic Quantization approach 2J to Quantum Field 
Theory and provide quite a general (and rather 
simple) way to perform Lattice Perturbation The- 
ory calculations in several contexts. In the work 
we are going to present in these pages, we ap- 
ply them to QCD at finite temperature, in the 
framework of a 3d effective theory. Lattice per- 
turbation theory calculations are needed in order 
to connect non-perturbative lattice results to the 
continuum. Hence we set up our NSPT simu- 
lations and compute the plaquette up to g 8 in 
the 3d pure gauge theory: we reach the infinite- 
volume value of each coefficient of the series by 
extrapolating the results we get for lattices of dif- 
ferent sizes, and in particular we find a logarith- 
mic divergence in the last coefficient. 

1. PHYSICAL FRAMEWORK 

The topic is hot QCD, and in particular the 
Free Energy Density, or the pressure of the quark- 
gluon plasma: it represents a good observable to 
study the deconfinement phase transition [3] . The 
aim is to fill the gap between perturbative meth- 



ods 0] and 4d finite-temperature lattice simula- 
tions 5 . The former have to stay in the ex- 
tremely high temperature sector because of the 
poor convergence of the series, while the high- 
est temperatures used in lattice simulations are 
about 4-^5 times the transition temperature, be- 
cause of computational resource limitations. 

Dimensional reduction is a way to access the in- 
termediate regions. One integrates out the 'hard' 
as well as the 'soft' modes to get an effective the- 
ory for the 'ultrasoft' modes [H]. The result is a 3d 
SU(3) pure gauge theory, whose contribution can 
then be estimated with lattice Monte Carlo tech- 
niques. To complete the matching of 3d lattice 
results to continuum ones [7], one needs lattice 
perturbation theory computations. 

2. THE METHOD 

2.1. The Algorithm 

According to the Stochastic Quantization ap- 
proach, we sample the phase-space of the field 
theory randomly, according to a stochastic equa- 
tion, the Langevin equation, 
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along a new, non-physical, stochastic time t. This 
sampling is such that the average over the noise 
i] leads to the Feynman-Gibbs ensemble we need 
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for functional integration: 



2.67 



V[<p}0[<f>]exp{-S[cf>}} .(2) 



Hence, perturbation theory is performed by using 
the gauge-field-version of the Langevin equation, 



dtU v = [-iVS[U v ]-ir)]U v , 



(3) 



and by replacing the gauge fields U v with their 
perturbative expansion g k U^ in the cou- 
pling constant g. Then we solve numerically the 
resulting system of equations via discretization of 
the stochastic time t — nr. 

2.2. The Code 

For this purpose we set up, from scratch, a 
C++ code (taking a hint from the code written 
for APE machines in their own language). The 
flexibility of this language makes it possible to 
set up quite a general environment of classes and 
methods to handle lattice-structures (both 'link- 
like' and 'site- like'). The code, moreover, is in- 
tended for PC-cluster usage, relying on MPI lan- 
guage for node-communications. In order to im- 
prove code-parallelization through spreading sub- 
lattices over the nodes, some tricks are used, like a 
wide utilization of pointers and a smart allocation 
of physical memory especially suited to enhancing 
communication rates. 

2.3. Extracting data 

In practice, we evolve our system according to 
a discrete-stochastic-time version of the Langevin 
equation, keeping all the orders (in g) up to the 
one we are interested in, and we do that with dif- 
ferent time-steps. Moreover, choosing the time 
steps small enough to fall in the linear region 
allows us to extrapolate in a simple way to the 
continuum-stochastic-time case. All this, then, is 
repeated for different lattice sizes in order to take 
the infinite- volume limit. The preliminary results 
we present here are the result of 50 days of run- 
ning on a cluster made of 10 bi-processor Athlon 
MP2100 PC's. 

3. RESULTS 

We computed the weak-coupling expansion for 
the pure gauge plaquette V = -^-Tr(nrjL/j). Our 
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Figure 1. Benchmark: lattice-size dependence of 
the coefficient 

expansions are written in powers of j3~ l — 
The first two coefficients of this expansion are 
known analytically for arbitrary dimensions 8 , 
so we can use them as a benchmark for our code. 
In particular, since for the leading coefficient also 
the finite-size corrections are known, we can as 
well check each single finite-size measurement. In 
Fig. ^ we show the coefficient versus the 
lattice-size. We tried to fit our numerical results 
(the points with their error-bars in the figure) 
with different forms for the finite-volume correc- 
tions, and we found that very clearly data pre- 
fer the inverse- volume correction (the dashed line 
in the figure). Moreover, both the coefficient of 
the volume-dependence and the infinite-volume 
extrapolation are in good agreement with the an- 
alytical results: 

(1 - V) A p-t (L) = 2.667(1) - 2.8(3)L~ 3 , (4) 

(I - -p) a g a ! ytic (L) = 2.66667 - 2.66667L" 3 . (5) 



Also for the f3~ 2 coefficient data prefer the 
inverse-volume as the leading term in the 'effec- 
tive' finite-size correction. In fact, for different 
power law fits, the infinite-volume extrapolation 
is quite stable and, again, in good agreement with 
the analytical result: 



{l-V)f-S{L = oo) 



1.95(1) , 



(6) 



(I - V) a JLi ytic (L = oo) = 1.94862. (7) 
The next coefficient, f3~ 3 , is the first original 
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Figure 2. The logarithmic divergence: lattice-size 
dependence of the (3~ A coefficient 

result: 

{l-V)p-3{L = oo) =6.7±0.2. (8) 

In |U] , the authors estimate this coefficient rescal- 
ing the analogous 4d value (which is known, |10p. 
finding (f — V)g-3 ~ 7.02, which is in quite rea- 
sonable agreement with the value we found. 

As for the (3~ A coefficient, we expect that the 
lattice reflects in the volume-dependence the loga- 
rithmic divergence one finds in continuum pertur- 
bation theory II . Indeed we found (see Fig. [21 
that our data-fit clearly improves if we add a log- 
arithmic term in the interpolation law. Writing 
the (3~ A term as 

(1 - P) p _< (L) = £ ln(i 3 ) + c + E fc c k L-\ (9) 
our fits result in the following estimates: 

£ data = q g ± q 4 _ c data ^ 2 3 ± 5 , (10) 

where the usual difficulties in fitting a logarithmic 
term leads to a quite big error for the £ coefficient. 

Furthermore, the coefficients of these two log- 
arithms (the lattice one, of the lattice-volume, 
and the continuum one, of the cut-off of the 
dimensional-reduced theory) must be the same. 
Hence the £ coefficient we extrapolate with our 
simulations can be seen as an additional (and, it 
turns out, positive) check for the reliability of our 
method. In fact, our estimate £ data is in good 
agreement with the continuum result found in 
[TT]: £ cont = 0.9765. 



As a next step, we can plug this analytical value 
£cont m |- our j n orc [ er to better estimate the 
constant coefficient Co, for which we then get the 
following result: 



data 
C 



log fixed 



25 ± 2 . 



(11) 



CONCLUSIONS 



We use NSPT methods to compute the plaque- 
tte in a 3d pure gauge SU(3) theory, up to g 8 . 

Up to now, we are still collecting statistics in 
order to improve our measurements, both by in- 
creasing the number of finite stochastic-time sim- 
ulations, and by filling the gaps over the range of 
lattice-sizes probed. 

Furthermore, our code is now ready to simulate 
an additional adjoint Higgs field coupled to the 
gauge sector, which comprises the effective theory 
for 'soft' modes in the plasma. We are performing 
preliminary simulations for measuring both the 
quadratic and the quartic Higgs condensates. 
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